Band engineering in graphene with superlattices of substitutional defects 
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We investigate graphene superlattices of nitrogen and boron substitutional defects and by using 
symmetry arguments and electronic structure calculations we show how such superlattices can be 
used to modify graphene band structure. Specifically, depending on the superlattice symmetry, the 
structures considered here can either preserve the Dirac cones (Dgh superlattices) or open a band 
gap (Z?3h). Relevant band parameters (carriers effective masses, group velocities and gaps, when 
present) are found to depend on the superlattice constant n as l/n p where p is in the range 1 — 2, 
depending on the case considered. Overall, the results presented here show how one can tune the 
graphene band structure to a great extent by modifying few superlattice parameters. 



I. INTRODUCTION 



Single layer graphene is a very promising material for 
a future silicon-free nanoelectronics. The peculiar char- 
acter of its charge carriers comes from the intersection of 
the 7r/7r* electronic bands occurring at the corners of its 
hexagonal Brillouin zone. This gives rise to the so-called 
Dirac cones at the Fermi level and makes graphene a 
zero-gap semiconductor in which low ener gy e xcitations 
behave as massless, chiral Dirac particles® 3 . In turn, 
this implies a series of interesting physical effects that 
open new perspectives for fabricating novel electronic 
devices^, e.g. high -per formance transistors for radiofre- 
quency applications^. In this perspective the possibility 
of engineering graphene's band structure by introducing 
defects, strains or external potentials has gained impor- 
tance in the recent past, in particular for opening a gap 
in the band structure which is essential to design logic 
devices. Indeed the non-vanishing residual conductance 
of intrinsic graphene avoi ds t he complete current pinch- 
off in the pristine material^, thereby limiting the on-off 
ratio to ~ 10 1 — 10 2 . A number of controlled techniques 
for energy band engineering have been proposed other 
than the actively pursued goal to obtain nanoribbons of 
controlled size and edge geometry. Most of them are 
based on the use of superl attice s of external potential d 9 * 10 ! 
or defects such as hore d 11 * 12 ^ and adsorbates^. Con- 
trolled vacancies on graphene^, as well as large holes 
symmetrically arranged to form graphene anti-dots^ 
have actually been realized with modern lithographic 
and self-assembling techniques. Preferential sticking of 
atom s indu ced by Moire patterns^ or by other electronic 
effect ^121121 could also induce a superlattice ordering that 
modifies graphene energy bands. Likewise, there is a 
great hope that novel bottom- up techniques^ may be ap- 
plied to fabricate atomically precise graphenic structures 
as already shown for nanoribbon& 2 °. These approaches 
might allow to realize in the near future graphene-related 
two-dimensional materials with modified characteristics, 
e.g. linearly-dispersing bands with variable Fermi ve- 
locities or semiconducting structures. In this paper we 
focus on atomically precise superlattices of substitutional 



atoms. The present work connects to and extend a recent 
worlJ^ where we have shown that in properly designed 
superlattices of holes or adatoms one can open a gap with- 
out breaking graphene point symmetry, i.e. preserving 
the pseudo-relativistic behaviour of charge carriers which 
makes graphene so attractive. The structures suggested 
in RefpD have n vacancies (hence missing p z orbitals) at 
the sites of a honeycomb superlattice, as a consequence 
of the introduction of C vacancies (holes) or chemisorp- 
tion of simple atomic species. Here we consider similar, 
highly symmetric structures, but with ir vacancies re- 
placed by boron and nitrogen. Similar defects have been 
recently considered for tuning the electronic properties of 
graphene nanoribbons and other carbon based structures 
suggesting that, when arranged to form particular struc- 
tures, they can turn the material into a semiconductor or 
a half-metaP^HH] Half-metallicity and the other many- 
body effects in such a structures open new perspectives 
in the field of carbon-based material s for spintronic ap- 
plications: for a recent review see Ref! 25 * 26 ! and references 
therein. 

In this paper we show that, depending on the super- 
lattice symmetry, one can obtain either electron (hole) 
doped substrates with pseudo-relativistic massless carri- 
ers or semiconducting structures with a quasi-conical dis- 
persion, and with the help of electronic structure calcu- 
lations (tight-binding and density-functional theory) we 
determine carriers velocities, effective masses and band 
gaps (when present) as functions of the superlattice . The 
focus is on boron and nitrogen, mainly because of the 
fast progresses in methods for the controlled synthesis of 
B and N doped graphemes. For instance, Panchakarla 
et alW^ have recently shown how it is possible to insert 
B or N dopants in graphene by adding the correct pre- 
cursors in the arc discharge chamber, while Ci et alW^ 
have reported the synthesis of large islands of boron ni- 
tride embedded in graphene by atomic layer deposition 
techniques. Methods to selectively replace C atoms in 
the graphene lattice have also been proposed^, thereby 
suggesting that the superlattice structures considered in 
this paper might soon become feasible. 
The paper is organized as follows. In the next section we 
summarize the computational details of the calculations. 



2 



Table I: Parameters used in the tight-binding Hamiltonian. 
All the values are in eV. 



Atom e 


ti 


t 2 


h 


C 0.000 


-2.900 


+0.175 


-0.155 


B -1.5225 


+1.450 






N +1.5225 


-1.450 







enough to reproduce the bands close to the Fermi en- 
ergy. 

Therefore, we computed DFT band structures for nxn 
graphenc superlattices up to n=14, and for larger struc- 
tures we relied on TB calculations only. 



III. RESULTS AND DISCUSSION 



Then, we show how p-(n-) doped graphene-like struc- 
tures result when substitutional defects are arranged in 
honeycomb superlattices, whereas semiconducting struc- 
tures with quasi-conical dispersion (massive Dirac carri- 
ers) result cither from a hexagonal supcrlatticc of from 
a honeycomb co-doped superlattice. Finally, we summa- 
rize and conclude. 

Throughout this paper we define the superlattice pe- 
riodicity using Wood's notation, i.e. by multiplying 
graphene's two-dimensional lattice vectors by the inte- 
ger (superlattice) constant n. 



II. COMPUTATIONAL METHODS 

The results shown in the next sections have been ob- 
tained from both tight-binding (TB) and density func- 
tional theory (DFT) electronic structure calculations. In 
the first case we diagonalized the usual tight-binding 
Hamiltonian for graphene it — tt* system, applying pe- 
riodic boundary conditions and including hopping terms 
up to the third nearest-neighbors. The on-site energies 
€i and hopping terms t\, <2 and £3 (for nearest, next-to- 
nearest and next-to-next-nearest neighbors, respectively) 
are those proposed by Nanda et aZPA They were fitted 
to accurate all-electron calculations to correctly repro- 
duce the Fermi velocity of single layer graphene. For 
the dopant atoms we only considered hoppings to near- 
est neighbor sites. Their values (ti), as well as those of 
the on-site energies (e»), are those introduced by Peres et 
aZpS, who have already successfully used them to study 
electronic effects in doped graphene. A summary of the 
TB parameters is listed in [IJ 

First principles DFT calcu lation s were preformed with 
the help of the VASP suit d 32 * 33 1 , using a supercell ap- 
proach. Core electrons were taken into account by projec- 
tor augmented wave (PAW) pseudo-potentials while for 
valence a 500 eV plane wave cutoff was used. To correctly 
represent the defect induced charge inhomogeneities 
we used the Perdew-Burke-Ernzerhof (PBE) gradient- 
dependent exchange and correlation functional 34 . Band 
structures were sampled by a T centered fc-points grid, 
never sparser than 6x6x1 in order to include every special 
point in the Brillouin zone (BZ). 

The TB parametrization was tested by comparing the 
band structure of few superlattices along the T-K-M- 
K'-r path with accurate DFT results. In every case 
the adopted parametrization was found to be accurate 



Graphene's peculiar electronic structure is strictly re- 
lated to the point symmetry of its lattice, D 6 h in the 
Shonflies notation. In the Brillouin zone, for each Bloch 
electronic state with k vector k, the relevant symmetry el- 
ements are those which either leaves k invariant or trans- 
forms it into one of its equivalent images, i.e. k — > k + G 
being G a reciprocal lattice vector. These elements form 
a subgroup of D^h, known as little co-group or simply 
fc-group at W 3 ^, which determines the possible symme- 
tries of the electronic states at k. At the high symmetry 
point K (or K') of graphene's Brillouin zone the fc-group 
is D^h, and Bloch functions built as linear combinations 
of p z orbitals span a two-dimensional irreducible repre- 
sentation (irrep) of such a symmetry group (E"). This is 
enough for the 7r — tt* degeneracy and the unusual linear 
dispersion at K (K 1 ). That this occurs exactly at the 
Fermi level is a consequence of the electron-hole (e — h) 
symmetry which approximately holds in graphene. In- 
deed, thanks to this extra symmetry, energy levels are 
always symmetrically arranged and, at half-filling, the 
Fermi level lies exactly at the center of the spectrum, 
where any doubly degenerate level is forced to lay^. In 
general, the number of doubly degenerate irreps in the 
BZ determines alone the presence of states (absence of a 
gap) at the Fermi level. We have recently showrPS how 
one can turn such number to be even at every special 
point -thus opening a gap in the band structure- by sym- 
metrically removing ll p z orbitals" in forming certain nxn 
superlattices. Substitutional defects behave similarly to 
p z vacancies (to which they reduce when the hoppings 
become zero) but introduce impurity bands which par- 
tially hybridize with those of the substrate. In addition, 
the diagonal disorder they introduce breaks e — h sym- 
metry giving rise to a Fermi level shift, i.e. to p— and 
n— doping for group IIIA and VA elements, respectively, 
as recently shown for both graphene^ and nanotubes^ 
. In the weakly defective superstructures considered in 
the following the defect-induced perturbation affects the 
electronic structure close to the Fermi level. With homo- 
geneous doping the latter shifts at most proportionally to 
1/n, i.e. as the square root of the defect concentration, 
as a consequence of the linear-energy dispersion which 
implies Ep = VFy/^n ei where Vf is the Fermi velocity 
of pristine graphene and n e is the electron (hole) excess 
density, n e oc 1/n 2 . Hence, analogously to the super- 
lattices of p z vacancies^, we make use of symmetry ar- 
guments to establish whether degeneracy occurs at the 
BZ special points in the important low-energy region. 
It is worth noticing at this point that, however small 



Figure 1: Folding of graphene Brillouin zone (BZ, blue line) 
into the superlattice ones (red filled hexagon) for some nxn 
structures, a) n=3ra, b) n—3m + 1 and c) n=3m+2, along 
with the case of y/3nxy/3nR30° superlattices (d). The K 
point in graphene's BZ, is labeled with a black filled dot. 



the defect perturbation is, the folding of graphene band 
structure occurs differently according whether the super- 
lattice constant n belongs to the sequence n = 3m or 
n = 3m + 1, 3m + 2 (m integer). As shown in [I] for the 
superlattices considered in this work (a-c), for n = 3to+1 
(3rn + 2) K and K' fold separately into K n (K' n ) and 
K' n (K n ), whereas for n = 3m they both fold to the BZ 
center r„. This means that n — 3m superlattices are 
expected to have rather unique properties related to the 
highly degenerate nature of the unperturbed spectrum. 
In the following we mainly focus on n = 3m + 1, 3m + 2 
superlattices and only occasionally look at the properties 
of n — 3m ones. A further six-fold superlattice symme- 
try, the v / 3^xv / 3ni?30° case reported in [ijd), will not 
be considered here since in that case band folding occurs 
analogously to the 3mx3m case discussed above. 



A. Honeycomb superlattices 

A honeycomb-shaped superlattice is a natural choice 
for nxn superlattices (nxn-honeycombs thereafter), since 
it preserves the D^h point group symmetry of pristine 
graphene. The superlattice unit cell contains two sub- 
stitutional atoms and is shown in [2j If the atomic radii 
of the dopants are small enough that lattice distortions 
are minimal, the system overall symmetry is preserved 
and Dirac cones at K n and K' n are expected. This is 
the case of boron and nitrogen substitutional defects, 
whose DFT-optimized structures show no appreciable 
lattice distortion. Both TB and DFT calculations con- 
firm that n = 3m + 1 and 3rn + 2 honeycomb superlat- 
tices made of B or N substitutional defects only show a 



Figure 2: A 4x4-honeycomb superlattice: the black line repre- 
sents the unit cell while the Wigner-Seitz and Brillouin zone 
are shown in yellow and green respectively. Red balls are 
sublattice substitutional defects forming the superlattice. 



low-energy band structure very similar to that of perfect 
graphene, but with the Fermi level lying respectively be- 
low (p-doped) and above (n-doped) their Dirac point. 
In principle, with properly designed n— or p— back- 
doping, e. g. ele ctric-field induced but also via molecular 
adsorptio n 1 38 1 39 1 , such shift can be offset and the analogy 
with pristine graphene can be fully exploited. 
[3]ja) shows the TB and first principles band structures 
of one nxn-honeycomb together with the position of the 
Fermi level ^c)) in such n- and p-doped superlattices at 
different impurity concentrations. As expected, the shift 
(A) of the Dirac cones with respect to the Fermi level 
(see[3jc)) is, to a good approximation, inversely propor- 
tional to the dopant concentration for both B and N dop- 
ing, though with opposite sign. The difference between 
TB and DFT band structure is minimal, and this con- 
firms that the tight binding parameters adopted are good 
enough for accurately describing the low-energy features 
of the nxn-honeycomb superlattices investigated in this 
paper. In[3]jb) we also report the unique band structure 
resulting from the special folding in the n = 3m sequence; 
as it is evident from the inset of ^b), the four- fold de- 
generacy occurring at r„ is partially lifted, and a gap is 
introduced in one of the two cone replica. 

The group velocity for electrons and holes taken close 
to the cone apex (but rather adequate for a wider en- 
ergy range) is shown in [3jcl) for n- and p-doped super- 
lattices. The two curves approach the limit of clean 
graphene with different trends. Upon non-linear curve 
fitting the group velocity v (relative to the one in pris- 
tine graphene) for p-doped honeycombs is found to be- 
have as v/vf oc 1 — n~ 1,29 , while for n-doped honeycombs 
as oc 1 — n -184 . The difference between the two cases is 
due to the value of the on-site energies and hopping of the 
dopants which determine the degree of hybridization of 
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Figure 3: a): TB (black lines) and DFT (red lines) band 
structures for the 4x4-honeycomb boron superlattice. b) the 
TB band structure of the 3x3-honeycomb boron superlattice 
arising from folding in n — 3m superlattices. The inset shows 
a close-up of the region close to T n . c) Absolute shift of 
the Dirac cones apex (A) with respect to the Fermi level, in 
p-doped (B, red) and, n-doped (N, black) honeycombs, d) 
Group velocity for charge carriers close to the cones apex 
for the boron (black) and nitrogen (blue) case. Circles and 
squares for n = 3m + 1 and n — 3m + 2. 

their impurity levels with that of bulk graphene. With 
the parameters used (see [I]) , which are symmetric with 
respect to the on-site energy of C atoms, this can only 
happen because of the asymmetry in graphene electronic 
structure introduced by the next-to-nearest neighbor in- 
teractions. Other superlattices made of group IIIA (Al, 
Ga, In) and VA (P, As, Sb) dopants have been tested 
by first principles calculations. In any case we found 
that, after geometric optimization of the lattice struc- 
ture, the impurities stand out from the graphene layer 
plane and considerably distort the neighboring lattice po- 
sitions. The resulting band structures are metallic but 
lack of Dirac cones due to the reduced symmetry. 



B. Hexagonal superlattices 

When one defect per supercell only is introduced a 
nxn hexagonal superlattice (a "nxn-hexagon" ) results, 
as shown in|4j This kind of structures is closely related 
to the honeycomb ones, having one extra substitutional 
atom at the center of a hexagon of defects. A closer in- 
spection, however, reveals that, due to the presence of the 
underlying C network, the point symmetry is reduced to 
-D3/1, with a planes missing with respect to the honey- 
comb counterparts. It follows that the fc-group at K n 
{K' n ) is C^h, with no irreducible two-dimensional (com- 
plex) representations (see km). Hence, degeneracy is re- 



Figure 4: A 2x2-hexagon superlattice: the black line rep- 
resents the unit cell boundary while the Wigner-Seitz and 
Brillouin zone are filled in yellow and green respectively. Red 
balls are substitutional defects positions. 
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Figure 5: Wigner-Seitz (yellow, left side) and Brillouin zone 
(right side) for the hexagonal superlattices. Considering the 
defects only the point group is Deh, and the corresponding 
fc-group in K n D^h (upper panel). Overall the underlying 
carbons remove the a 1 planes, reducing the symmetry to Dzu- 



moved at the special points and a (small) gap opens in the 
band structures, close to the (shifted) Fermi energy. This 
is shown in[6ja) where the TB and DFT band-structures 
of the 4x4 hexagon are reported. The energy spectrum of 
such gapped graphene is compatible with charge carriers 
behaving as massive Dirac particles 

E(k) = ±v^k% + k 2 y + m* 2 v 2 (1) 




Figure 6: a) TB (black) and DFT (red) band structures 
for the 4x4-hexagon. b) the TB band structure of the 3x3- 
hexagon boron superlattice arising from folding in n = 3m 
superlattices. The inset shows a close-up of the region close 
to r„. Carriers effective speeds (c), masses (d) and energy 
band-gaps (e) versus 1/n. Black and blue values refer to B 
and N superlattices, respectively. 

where v,m* are the effective 'speed of light' and 'rest 
mass', respectively, and determine the gap size 

AE = 2m* 'v 2 (2) 



Figure 7: A 4x4-honeycomb superlattice: the black line repre- 
sents the unit cell while the Wigner-Seitz and Brillouin zone 
are shown in yellow and green respectively. Red and Blue 
balls are B and N substitutional defects. 

the massive, pseudo-relativistic energy dispersion is su- 
perimposed with a massless one, thereby giving rise to a 
two-valley system with very different charge carriers. As 
shown in the next subsection, all the features discussed 
in this Section can be brought at the Fermi level by co- 
doping the substrate in forming a honeycomb structure 
with the same D 3 h symmetry discussed here. 

C. Co-doped superlattices 



According to the semiclassical theory of conduction m* 
is also the effective mass rn e ff governing charge carrier 
mobility for k << mv\ for k » mv carriers behave 
pseudorelativistically with m e ff = and limiting speed 
v. The values v, m* and AE have been obtained by 
non- linear curve fitting of the numerical results to [T] and 
are reported in panels (b)-(d) of 6j For rtxn-hexagons the 
band gap is very dependent on the type of dopant ([6j(e)): 
the maximum gaps, occurring in 2x2 hexagons, are 0.93 
eV for nitrogen and only 0.17 eV for boron. The effec- 
tive masses of electrons and holes (|6](d)) roughly scale 
as the gaps: oc n~ 145 and oc n~ 152 for n- and p-doped 
structures respectively, and their maximum is 3.7xl0 -2 
and 6.7xl0 -3 m e . This is similar to the case of graphene 
nanoribborP&l, whose band gap scales as the inverse of 
their width even though here the gap is due to symmetry 
breaking rather than quantum confinement. The shift of 
the Fermi level (not shown) is again proportional to the 
square root of the defects concentration, that is now only 
half of the value for honeycombs with the same super- 
lattice periodicity. Charge carriers velocities scale simi- 
larly for the two dopant species as shown in [6j (c) with 
a best-fit exponent close to -2 (v(B)/vf oc 1 — rt -198 , 



}(N)/v F oc 1 - 



-2.28 



). In H3l panel (b), we also re- 



port the particular band structure arising in n = 3m 
hexagon superlattices. At the relevant special point r„, 



One further possible superlattice arrangement is ob- 
tained by using two different dopants in the nxn- 
honeycomb unit cell, i.e. co-doping the structures with 
boron and nitrogen (see FjPb. In this way B and N 
atoms form a boron nitride-like honeycomb superlattice 
in which sublattices equivalence (and symmetry) is bro- 
ken. This is analogous to place graphene in the mod- 
ulating field of a proper substrate, e.g a hexagonal BN 
(0001) surface, which has been shown to lift the degen- 
eracy of the 7r — 7T* bands^ similarly for deposition, or 
growth, on silicon carbide surfaceiP^ES. The superlattice 
structures considered here offer the possibility to mod- 
ify the periodicity of the perturbation, and thus to tune 
the gap. Indeed, this kind of superlattices present D^h 
point symmetry, hence a Czk fc-group in K n ,K' n , and, 
analogously to the hexagonal case discussed above, open 
a band gap typical of massive Dirac particles. Differently 
from before, however, the structures considered here are 
iso-electronic with graphene and therefore the gap lies 
exactly at the Fermi energy. [8] shows the computed band 
structure (panel (a)), together with the values of the ef- 
fective speed of light (c), effective mass (d) and band 
gaps (e), obtained as in previous section by fitting of 
the numerical results, for different BN nxn-honeycombs. 
The results confirm the expectations, and show that such 
structures present a band-gap at the Fermi energy, com- 
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Figure 8: a) TB and DFT band structures for the co-doped 
4x4-honeycomb. b) the case of 3x3 superlattice with the low- 
energy region in the inset, c) Effective speed, d) mass at 
the K n point versus 1/n, and e) band-gap. The red and the 
blue dashed lines are the values for n- and p-doped hexagons, 
shown for comparison. 



patible with pseudo-relativistic massive carriers. Their 
effective rest mass is rather small, scales as oc n~ 1Ae 
and it is never larger than 0.016 m e for n > 4. This 
value compares favorably with the effective masses in 
Bii_j;Sb K topological insulators (m*=0.009 m e f^, and is 
generally lower than in bilayer graphene (to* =0.03 m e J^ 
or in any other traditional bulk semiconductors, such as 
InSb (m* =0.016 m e ). Since to* is the main factor affect- 
ing carrier mobility, the suggested structures turn out 
to be a good compromise between the need of opening 
a gap for logic applications and the desire of preserving 
the high mobility of charge carriers. In |8j panel (b) we 
also report the band structure of the n = 3to case. The 
structure is that of a zero-gap semiconductor, with two 
distinct charge carriers: one of them behave as electron 
(hole) in graphene, showing typical effects expected for 
massless carriers; the other is a more conventional one, 
with a finite excitation energy across a gap. 



or the honeycomb one is symmetrically co-doped, the 
Dirac cones detach from each other to form a gapped, 
quasi-conical structure whose excitations correspond to 
massive Dirac particles. Note that this situation clearly 
differs form the case of randomly arranged B or N im- 
purities, in which the density of states shows no band 
gaiP. 

For zero-gap structures the use of this superlattices offers 
the possibility to control the Fermi velocity by changing 
the structure periodicity, thereby offering the opportu- 
nity to investigate its role in the charge transport prop- 
erties. Differently from our recent proposal, the gapped 
band structures arise because of symmetry breaking, as 
in the case of graphene interacting with a substrate such 
as SiC or BN. In the same fashion the band gap size de- 
pends on the superlattice periodicity. In our calculations 
we have found that gaps and charge carriers velocities 
effective masses depends on l/n p , where p is in the range 
1 — 2, hence on some small power 0.5 — 1 of the dopant 
concentration, and on the dopant type (B or N). Over- 
all the structures proposed here show a band gap larger 
than ksT at room temperature, with an effective mass 
generally lower that 0.01 m e for reasonably dense meshes 
(n=4-10). Thus, the new class of graphene structures 
proposed might be promising candidates for the fabri- 
cation of high performance interconnects, valley-based 
devices^, but also for logic transistors, where a band 
gap is needed, but the extraordinary properties of pris- 
tine graphene need to be preserved. 
The electronic properties of these impurities superlattices 
rely on symmetry, hence are necessarily sensitive to the 
dopant positions. As a consequence, an accurate control 
of the system geometry is necessary to exploit their prop- 
erties. This might be possible in the near future with pre- 
cise bottom- up techniques, such as the ones recently used 
by Ruffieux and co-worker d 19 * 20 ! to fabricate nanoribbons 
of well-defined widths and edges. An indication of the 
possible synthetic routes together with formation ener- 
gies for such defects superlattices can be found in the 
Supporting Information. This information is available 
free of charge via the Internet at http://pubs.acs.org] 



SUPPORTING INFORMATION 



IV. SUMMARY AND CONCLUSIONS 

To summarize, we have studied the effects substitu- 
tional defects such B and N species have on graphene 
electronic structure when they are periodically arranged 
to form some superlattices. Using group theoretical argu- 
ments and both TB and DFT calculations we have shown 
that defects can either preserve the Dirac cones or open a 
band gap, depending on the superlattice symmetry (D^h 
and Dsh, respectively). Specifically, honeycomb-shaped 
superlattices of B or N atoms give rise to p- and n-doped 
graphene, respectively, preserving the Dirac cones. On 
the other hand, when a hexagonal superlattice is formed, 



In order to compare relative structural stabilities, for- 
mation energies for superlattice structures have been 
computed as follows: 

AH form = E S l - (2n 2 - T))fx c - n N ^ N - ubUb (3) 

where Esl is the total energy of the structure, n is the 
superlattice constant, r\ is the number of dopants per cell 
(1 for hexagons and 2 for honeycombs), n.j and /i.; are the 
number of atoms and the chemical potentials for each 
species. The chemical potentials were computed with re- 
spect to single layer graphene, gaseous N 2 and cv-Boron 
(in the so-called R12 structure). Density functional the- 
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Figure 9: DFT formation enthalpies (in eV) for Boron 
(black), Nitrogen (red) and co-doped (blue) honeycombs. 
Dashed lines show non-linear fits as a guide for the eye 

Table II: DFT formation enthalpies for doped superlattices 
for structures up to n—8 arranged by doping and by super- 
lattice symmetry. All the values shown are in eV 



Hexagon Honeycomb 



n 


B 


N 


B 


N 


BN 


2 


1.840 


1.078 


3.780 


1.754 0.841 


4 


1.513 


0.878 


3.455 


1.871 


1.500 


5 


1.499 


0.933 


3.224 


1.622 


1.593 


7 


1.380 


0.676 


3.063 


1.596 


1.626 


8 


1.394 0.695 


3.005 


1.490 


2.023 



ory results are shown in |TT] and in figure [9] 



For all the cases considered the formation energy of the 
superlattice structures is endothermic with respect to the 
pure elements with opposite trends for n- (p-)doped and 
co-doped structures. While two impurities of the same 
kind gain energy when lying further apart to each other, 
B and N tend instead to cluster together, in accordance 
with recent experimental observations^. However, in or- 
der to cluster the impurities have to diffuse through the 
graphene lattice by exchanging its position with a carbon 
atom. Such a process has a large energetic barrier being 
the two atoms covalently bound, hence the kinetics of this 
process is expected to be extremely slow. We therefore 
suggest that these structures are expected to be stable 
at room temperature as well as other known structures 
(e.g. graphene itself with respect to diamond). 
About the practical feasibility of substitutional defects 
superlattices we expect they can be produced by the same 
bottom-up approach recently used t o fabr icate graphene 
structures with atomic-scale controt^oi. According to 
this method an appropriate polyphenylene precursor, 
e.g. produced by Ullmann coupling, undergo cyclode- 
hydrogenation on a Cu surface to form a polycyclic aro- 
matic hydrocarbon (PAHjSlHIID that in our case should 
include the dopant atom. When appropriately func- 
tionalized the PAHs can then polymerize on the metal 
surface to form graphene domains up to a nanometer 
scale as recently shown for the synthesis of the atomi- 
cally precise nanoribbons^l. Such a technique has been 
also adopted to produce nitrogen doped fullerenes^ and 
two-dimensional polymers^, hence it might indeed lead 
to graphene-based superlattices of substitutional defects 
starting from chemically doped polycyclic aromatic hy- 
drocarbons. 



* Electronic address: simone.casolo@unimi.it , Electronic ad- 
dress: rocco.martinazzo@unimi.it 
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